Pinwheel stability, pattern selection and the geometry of visual space 
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It has been proposed that the dynamical stability of topological defects in the visual cortex reflects 
the Euclidean symmetry of the visual world. We analyze defect stability and pattern selection in a 
generalized Swift-Hohenberg model of visual cortical development symmetric under the Euclidean 
group E(2). Euclidean symmetry strongly influences the geometry and multistability of model 
solutions but does not directly impact on defect stability. 



oo 
o 
o 

a 

<N 



d 

• i-H 

<N 
> 

CO 
00 

cn 
o 

00 

o 



X 



When an object in the visual world is rotated about the 
axis of sight the orientations of distant contours change in 
a coordinated manner, leaving relative orientations of dis- 
tant contours invariant. In the visual cortex of the brain, 
it has been predicted that this geometrical property of 
visual space imposes the so called shift-twist symmetry 
on joint representations of contour position and contour 
orientation [I]. This symmetry requires the equivariance 
of dynamical models for these representations under the 
Euclidean group E(2). In particular, it has been hy- 
pothesized that shift-twist symmetry stabilizes topologi- 
cal pinwheel-defects in models for the emergence of ori- 
entation selectivity during cortical development [2j[3j|4]. 
Pinwheel defects are singular points in the visual cor- 
tex around which each orientation is represented exactly 
once [5j [6j [7]. They are initially generated in the visual 
cortex at the time of eye opening |8, 9|. This process has 
been theoretically explained by spontaneous symmetry 
breaking [TQ| ITT] . Why in the brain, pinwheels remain 
present at all developmental stages, although they are 
dynamically unstable in many models of visual cortical 
development [TTJ [I2j HSl [14] and in analogous phys- 
ical systems [TSJ [161 E] •> remains unclear. 

Previous studies [£TJ [H E9 CE3 found pin- 
wheels dynamically unstable only in models exhibiting an 
E(2)xU(l) symmetry, which is higher than the E(2) sym- 
metry of visual perceptual space. In contradistinction, 
models exhibiting only Euclidean symmetry have been 
shown to exhibit stable pinwheels [2] [3j [18] suggesting 
that the stabilization of pinwheel defects may be closely 
related to the 'reduced' Euclidean symmetry. There are, 
however, also other scenarios that predict the emergence 
of stable pinwheels with higher than Euclidean symmetry 
[12] [14] leaving it hard to judge to actual role of Euclidean 
symmetry in orientation map development. 

Here we analyze the impact of Euclidean symmetry on 
pattern selection, i.e. the question of whether stable pin- 
wheel arrangement exist and what their geometric orga- 
nization is. We construct a generalized Swift-Hohenberg 
model [T9| [20] symmetric under the Euclidean group 
E(2) that allows to study the transition from higher 
E(2)xU(l) to lower E(2) symmetry by changing a param- 
eter that controls the strength of shift symmetry breaking 
(SSB). Using weakly nonlinear analysis we derive ampli- 



[e|=0 



< |e| < 1/2 



1/2 < \e\ 





w 





Figure 1: (a) Plane waves with wavevector in horizontal (top) 
and oblique (bottom) direction for variable strength of SSB 
(e = 0, 0.35, 1) (b) PWCs of varying intersection angle a, 
tt/4 < a < 7r/2. (c) Energy of solutions depends on e and 
on a. For sufficiently large e PWCs are energetically favored 
relative to plane waves, (dashed: energy of plane waves, plain: 
energies of PWCs for a = 7r/4, 7r/3, tt/2) 



tude equations for stationary planforms and find three 
classes of stationary solutions: stripe patterns without 
any pinwheels, pinwheel crystals with pinwheels regularly 
arranged on a rhombic lattice, and quasi-periodic pat- 
terns containing a large number of irregularly spaced pin- 
wheels. We calculate the phase diagram of these solutions 
depending on the strength of SSB, the effective strength 
of nonlocal interactions, and the range of nonlocal inter- 
actions. With increasing strength of SSB, pinwheel free 
patterns are progressively replaced by pinwheel crystals 
in the phase diagram while both pinwheel free patterns 
and pinwheel crystals remain stable. Phases of aperiodic 
pinwheel rich patterns remain basically unaffected. A 
critical strength of SSB exists above which multistable 
aperiodic patterns collapse into a single aperiodic state. 
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The spatial structure of an OPM can be represented 
by a complex field z(x) where x denotes the 2D position 
of neurons in the visual cortex, 0(x) = arg (z(x)) /2 their 
preferred stimulus orientation, and the modulus \z(x)\ is 
a measure of their selectivity [10]. In this representation, 
pinwheel centers are zeros of the field z(x). The sim- 
plest models for the formation of OPMs are defined by a 
dynamics 



dtz(x,t) = F[z](x,t). 



(1) 



where t denotes time and F[z] is a nonlinear operator. 
We assume the dynamics equivariant under translation 
T y z(x.) = z(x + y), rotation R a z(x.) = e 2ta z(Q- a x.) 
with rotation matrix f^ 5 and reflection at the corti- 
cal (1,0) axis Pz{x) = z(~x), thus expressing the fact 
that within cortical layers there are no special locations 
or directions |21j. In addition, if interactions between 
OPM development and visuotopy are neglected it is also 
equivariant under global shifts of orientation preference 
Spz(x) = e z/3 z(x) (shift symmetry). Rotations R thus 
consist of a composition of phase shifts S and coordinate 
rotations D, i.e. R a — S2 a °D a with D a z(x) = z(Q- a x.). 

We consider the general class of variational models |3, 
[TD1H3] for which F[z] has the form 



F[z] = Lz + eMz + N 3 [z}. 



(2) 



Here L is a linear, translation invariant and self-adjoint 
operator, that accounts for a finite wavelength instability. 
7V3 is a cubic nonlinearity which stabilizes the dynamics. 
The second term involves a complex conjugation Cz — 
z and thus manifestly breaks shift symmetry when e 7^ 
0. M is assumed to be linear, translation invariant and 
bounded. Equivariance under rotations, [MC, R a ] = 0, 
requires 



and equivariance under parity [M, P] = 0. 

As a concrete example we will consider the model 
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L = r -(k z c +V z ) 
N 3 [z] = (l- 5 )|z(x)| 2 z(x) 
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(3) 

(4) 
(5) 



x J d 2 y (\z(y)\ 2 z(x) + \z{y) 2 z{^j ^ 
M = r{d x +id y )\d xx +d yy )- 2 (6) 

where L is the Swift-Hohenberg operator [19j [20] with 
critical wavenumber k c and instability parameter r. N3 
is adopted from [14], where a sets the range of the non- 
local interactions and g determines whether the local 
(g > 1) or the nonlocal term (g < 1) stabilizes the dy- 
namics. M is the simplest differential operator which 
transforms according to Eq.(|3|. It is unitary with spec- 
trum oc e 4iarg ( k ). 

We used weakly nonlinear analysis [22] to study poten- 
tial solutions of Eq.([2|. We consider planforms z(x) = 



Sj=o 1 ^je* kj ' x , kj = /c c (cosaj, sinaj) with 2n modes 
where we require that to each mode also its antiparallel 
mode is in the set. By symmetry, the dynamics of the 
amplitudes Aj at threshold has the form 



2n-l 



2n-l 



Aj — A.j~\~€.A.j 6 
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^ gjk\A k \ 2 Aj- ^ f jk A k A k _Aj 



k=0 



(7) 

where j denotes the index of the mode antiparallel to 
mode j, kj- = — kj and with real valued and symmet- 
ric matrices gj k and fj k which determine the coupling 
and competition between modes. They can be expressed 
in terms of angle-dependent interaction functions g(a) 
and f(a), which are obtained from the nonlinearity Ns[z] 
(cf.p~4j HU [22]) • For simplicity we restrict the following 
analysis to the class of permutation symmetric models, 
defined in [14], for which g(a) = g(a + it). [23] 

The simplest solution to Eq.([7f is obtained for n = 
1 and consists of plane waves with wavevector k = 
/c c (cosa, sin a). For |e| < 1/2 it is given by 



z(x) = — — [y/l + 2e cos(kx + 0)+iV / l - 2esin(kx + (/>)] 
V^oo 

(8) 

with arbitrary phase <j). Hence with SSB orientation an- 
gles are no longer equally represented. For e > cortical 
area for orientations a and a + 7r/2 is recruited at the ex- 
pense of a + 7r/4 and a + 37r/4 (and vice versa for e < 0). 
Beyond a critical strength of SSB, e* = 1/2, patterns 
only contain two orientations, z(x) = e 2%(x M cos(kx + <p) 
for e > e* and z(x) = ie 2%cx M sin(kx + <p) for e < — e* 

with M =V4(l+Vl)/30oo (Fig|J). 

Another class of solutions, rhombic pinwheel crystals 
(PWCs), exist for n = 2 and consist of two pairs of 
antiparallel modes forming an angle < a < tt/2 which 
are characterized by \Aq\ = \A -\ = a = \Ai\ = |^4i-|. 
We consider w.l.o.g. the case ao = —a/2 and a\ = a/2 
(FigjlJ)). With A 0y i = ae^ ' 1 , A - A - = ae iVQ ^ and 
Eo,i : = Mo,i + ^0,1 the stationary state is given by 
Ei' = -E and a 2 = (1 + ecos(E + 2a)) /C where C = 
3^oo + 2^oi+2/oi cos2Eo. The phase Eo is the solution to 
= sin 2E + e[sin(E - 2a) - (2 + 3g 00 /goi ) sin(E + 2a)] 
which bifurcates from Eo = ±7r/2 for e = 0. The energy 
is given by E PWC = -4a|(l + |e| cos(E + 2a)) + 2a 4 (. 
For the model Eqs. ppl ) the e and a dependence of the 
energy is shown in Figjlp. The solution then reads z(x) = 
2a [e iS °/ 2 cos(k x + A /2) + e'^' 2 sin(kix + Ai/2)] 
with arbitrary Ao and Ai. 

A large set of quasiperiodic solutions originates 
from the essentially complex planforms (ECP) z(x) = 



^i=c ) A j e that solve E( l- Y] for e = O* Here ' 
wave vectors kj = k c (cos j^j, sin^j) (j = 0, . . . , n — 1) 

are distributed equidistantly on the upper half of the 
critical circle and binary variables lj = ±1 determine 
whether the mode with wave vector kj or with wave vec- 
tor — kj is active (FigMa) left column). 

For |e| > we findthat ECPs generalize to z(x) = 
Y^JjZl [A~j~e tlj ^ jyi + Aje~ ll ^ j *] . For each n there exists 
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Figure 2: Attractors of amplitude equations (Eqj7| with full, 
partially broken and completely broken shift symmetry, (a) 
ECPs. Preferred orientations are color coded [see bars in 
(b)]. Arrangement of active modes on the critical circle and 
corresponding OPMs. For n = 3 and 8 there are 2 and 15 dif- 
ferent classes of ECPs, respectively. Complete (partial, no) 
suppression of opposite modes for full (weakly broken, max- 
imally broken) shift symmetry (left, middle, right column). 
(b) OPM in tree shrew VI (data: L.E.White, Duke Univ., 
USA). Arrows pinwheel centers. Scale bar 1 mm. (c) With 
increasing degree of symmetry breaking e amplitudes of an- 
tiparallel modes A- grow and eventually (at e = e*) reach the 
same absolute value as active modes A+. 
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: 7 [#oo + 2 Ej=o fij cos where 
When |e| < |ej stationary ampli- 



tudes 
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fulfill 



7 (1 ± y/l-e 2 /el). When 



|e| > |ej the amplitude of antiparallel and active modes 
are equal a\ = j[l + |e|]/[l + e*] (Figj^j:). A simple 
measure of the degree to which SSB affect n-ECPs is 
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For a stationary n-ECP we find q 
and q = signe if |e| > |ej. 



Stationary phases (jr^ 
i>j + <\>~ = (47r/n) j if e > and <j>f + <j>J = (4?r/n) j + 
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Figure 3: (a)-(c) Pinwheel densities for all realizations of 
ECPs with 3 < n < 17 and different degrees of shift symmetry 
breaking e. (d) Pinwheel densities for n — 17 (dots) and for 
n — > oo in the Gaussian approximation (gray region). 



if e < 0. This implies that all orientations are represented 
in patterns with n > 3. The solution can be written 

n-l 

*(x) = v^^Iv/TT^x,^) + y/T^zfaM 



3=0 



with 



cos(/jkjX + <j>j) and 2?(x, <^-) 



ie z ~^ J sin(/jkjX + 0j) and arbitrary phases <pj. For n = 1 
this is in agreement with Eq.(|8| where e* = 1/2. Re- 
flection at the axis parallel to kj acts on the functions z e - 
and z° as +1 and —1, respectively. Thus z e - and z° corre- 
spond to the even and odd eigenvectors of the nullspace 
of L + eMC (cf.HBl). For e > (e < 0) the even (odd) 
part dominates the solution. 

The dynamics Eq.([7]) exhibits a potentially exceedingly 
high number of multistable solutions. The energy of n- 
ECPs is given by E n = — 727 [1 + e 2 /e*] for |e| < e* and 
E n = —717(1 + |e|) 2 /[l + e*] for |e| > e*, respectively, and 
does not depend on the variables lj which identify a par- 
ticular n-ECP. Due to the growth of antiparallel modes 
with increasing |e| patterns for all different realizations 
lj with phases 4 



l 3®3 



i(l-sign(e))(l-Z> 
arbitrary but fixed) eventually collapse in a single state 

*(x)«E£o*i /o ( x >&) ( F W- 

This collapse manifests itself in the pinwheel densi- 
ties p n of various ECPs shown in Fig.3a-c for different 
strength of SSB. The pinwheel density of an n-ECP in 

and 



7T 



the large n limit is p(x) — 1 — ^-(1 — e 2 /e 2 )(" 2 

depends through £ = ^|Ej=o 'jkjl ^ 1 011 the con- 
figuration of active modes. Fig. 3d shows that pinwheel 
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Figure 4: Phase diagrams of the model, Eqs.(2]|6), near crit- 
icality for variable SSB e. The graph shows the regions of 
the g — a/A plane in which n-ECPs and PWCs have minimal 
energy (n = 1 — 25, n > 25 dots). Regions of maximally bro- 
ken shift symmetry [e > e* (iV, g, a)] shaded in gray. Regions 
where a-PWCs prevail is shaded in blue, intensity level codes 
for the relative angle a. (light blue: tt/4 < a < tt/2 -.dark 



densities fill a band of values. With increasing degree of 
SSB this band narrows and pinwheel densities eventually 
equal 7r at the critical value |e| = e*. 

To reveal how SSB affects pattern selection we calcu- 
lated the phase diagram for the model specified in Eqs.(|4| 
[6| for various values of e. FigJI] shows the configura- 
tions of n-ECPs and PWCs minimizing the energy |24j . 
Plane waves are progressively replaced by a-PWCs with 
increasing SSB. Depending on the location in parameter 
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Figure 5: Stability regions of ECPs with n — 20 for different 
e. Dashed line denotes the critical line e = e*(iV, cr), above 
which q = 1. Shaded region denotes region in the g — a/A 
plane for which that planform is a stable solution of the dy- 
namics and coexists with planforms of nearby values of n, e.g. 
n = 18, 19, 21, 22. In the inner region (marked by the two 
inner lines) this solution minimizes energy. 



space and on e, a particular angle a minimizing the en- 
ergy (c.f. Fig|lj;). Large n— ECPs are selected when the 
dynamics is stabilized by long-range interactions (g < 1, 
a > A). In this parameter regime plane waves and pin- 
wheel crystals are unstable. The degree of SSB q manifest 
in a given n— ECP attractor depends on e and on the lo- 
cation in the phase diagram. Above a critical line defined 
by |e*(iV, g, <j)\ = \e\ antiparallel modes are maximal and 
\q\ = 1 (gray area), below that line \q\ < 1. Figs. [2] and [5] 
show the high sensitivity of the dynamics to even small 
amounts of SSB, a substantial area in phase is occupied 
by ECPs with \q\ = 1 even for e = 0.02. 

Our analysis of pattern selection in visual cortical de- 
velopment demonstrates that dynamical models of orien- 
tation map development are very sensitive to the presence 
of interactions imposed by Euclidean E(2) symmetry. It 
also reveals that the impact of the Euclidean symmetry of 
perceptual visual space is not an all or none phenomenon. 
For weak SSB our Euclidean model closely mimics the be- 
havior of models possessing the higher E(2)xU(l) sym- 
metry. The only qualitative change that we found in the 
transition from higher to lower symmetry was the col- 
lapse of multistable solutions. This collapse, however, 
only happens when a finite critical strength of SSB is 
reached. Up to this threshold strength, models exhibiting 
E(2) and E(2)xU(l) symmetry seem to be topologically 
conjugate to one another. 

Our analysis also reveals that the impact of Euclidean 
symmetry differs qualitatively for aperiodic and peri- 
odic patterns. For aperiodic patterns, (i) the parameter 
regime in which they possess minimal energy is virtually 
unaffected by the strength of SSB, (ii) sets of different 
multistable solutions become progressively more similar 
and finally merge forming a single unique ground state, 
and (iii) Euclidean symmetry geometrically manifests it- 
self in specific two point correlations that are absent in 
the E(2)xU(l) symmetric limit. For periodic patterns 
such as pinwheel crystals and pinwheel free states, (i) 
the parameter regime in which these solutions possess 
minimal energy depends on the strength of SSB, and (ii) 
Euclidean symmetry geometrically manifests itself in a 
specifically selected tilt angle of rhombic pinwheel crys- 
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tals and a wave vector dependent underrepresentation of 
particular orientations for pinwheel free states. 

The Swift-Hohenberg model of Euclidean symmetry 
considered here predicts that aperiodic pinwheel rich pat- 
terns resembling the architecture of the primary visual 
cortex are only stable when long-range interactions dom- 
inate pattern selection, confirming previous predictions 



of a model of higher E(2)xU(l) symmetry. Our analy- 
sis predicts that in this regime, two point correlations 
provide a sensitive measure of the strength of SSB. Such 
correlations are therefore a promising tool for probing the 
impact of Euclidean symmetry in visual cortex develop- 
ment. 
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